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Isotropic-genesis nematic elastomers (IGNEs) are liquid crystalline polymers (LCPs) that have 
been randomly, permanently cross-linked in the high-temperature state so as to form an equilib- 
rium random solid. Thus, instead of being free to diffuse throughout the entire volume, as they 
would be in the liquid state, the constituent LCPs in an IGNE are mobile only over a finite, seg- 
ment specific, length-scale controlled by the density of cross-links. We address the effects that such 
network-induced localization have on the liquid-crystalline characteristics of an IGNE, as probed 
via measurements made at high temperatures. In contrast with the case of uncross-linked LCPs, 
for IGNEs these characteristics are determined not only by thermal fluctuations but also by the 
quenched disorder associated with the cross-link constraints. To study IGNEs, we consider a micro- 
scopic model of dimer nematogens in which the dimers interact via orientation-dependent excluded 
volume forces. The dimers are, furthermore, randomly, permanently cross-linked via short Hookean 
springs, the statistics of which we model by means of a Deam-Edwards type of distribution. We show 
that at length-scales larger than the size of the nematogens this approach leads to a recently pro- 
posed, phenomenological Landau theory of IGNEs [Lu et al., Phys. Rev. Lett. 108, 257803 (2012)], 
and hence predicts a regime of short-ranged oscillatory spatial correlations in the nematic align- 
ment, of both thermal and glassy types. In addition, we consider two alternative microscopic models 
of IGNEs: (i) a wormlike chain model of IGNEs that are formed via the cross-linking of side-chain 
LCPs; and (ii) a jointed chain model of IGNEs that are formed via the cross-linking of main-chain 
LCPs. At large length-scales, both of these models give rise to liquid-crystalline characteristics that 
are qualitatively in line with those predicted on the basis of the dimer-and-springs model, reflecting 
the fact that the three models inhabit a common universality class. 

PACS numbers: 61.30.Vx,61.30.-v,61.43.-j 



I. INTRODUCTION 



Nematic elastomeric materials 



Nematic elastomers are fascinating materials in which 
the liquid crystalline order is strongly coupled to the elas- 
ticity of the underlying elastomeric network (see, e.g., 
Refs. [1-6]). This strong nemato-elastic coupling gives 
rise to novel, emergent properties in nematic elastomers 
that are found neither in liquid crystal nematics nor or- 
dinary rubbery materials. One well-known example of 
such properties is the soft elasticity characteristic of mon- 
odomain nematic elastomers (see, e.g., Ref. [7]) that are 
formed via the Kiipfer-Finkelmann procedure [8] . 

Not only are nematic elastomers fascinating but also 
they have proven to be a challenging subject for the- 
oretical investigation. Part of the challenge originates 
in the dependence of the physical characteristics of ne- 
matic elastomers on the conditions under which the elas- 
tomers are prepared. For example, isotropic-genesis ne- 
matic elastomers (or IGNEs) — nematic elastomers cross- 
linked in the isotropic state — exhibit the so-called super- 
soft version of elastic response at sufficiently low tem- 
peratures (see, e.g., Refs. [9-11]), unlike their nematic- 
genesis nematic elastomer (or NGNE) counterparts. Fur- 
thermore, in thermal equilibrium the nematic alignment 



in IGNEs exhibits a polydomain structure (see, e.g., 
Refs. [3, 12-16, 18]) characterized by a length-scale deter- 
mined solely via thermodynamic quantities such as the 
density of cross-links. 

Another challenge to theorists comes from the fact 
that nematic elastomers possess a multi-level hierarchy 
of interdependent elements of randomness. First, there 
is quenched disorder in its conventional form, associated 
with the permanent chemical structure that originates 
in the cross-linking process. Second, as a result of suf- 
ficient cross-linking there arise the mean positions and 
r.m.s. displacements of the spatially localized polymers 
that constitute the elastomeric network, both of these 
elements being random. Third, there is the thermal dis- 
order associated with the Brownian motion of the posi- 
tional and orientational (i.e., nematogenic) freedoms in 
the state of the system just prior to the instant of cross- 
linking; in part, this thermal disorder is frozen in via 
the process of cross-linking. Fourth, there is also the 
thermal disorder associated with the Brownian motion 
of the nematogens in the state of the system long after 
the instant of cross-linking. It is not a priori evident 
how the interplay between the various types of random- 
ness present in nematic elastomers resolve themselves, 
and thus determine the equilibrium structure and elastic 
response of such media. As we shall try to make clear 



in subsequent sections of the present Paper, in order to 
understand nematic elastomers it is valuable to go be- 
yond the conventional notion of quenched disorder and, 
instead, to consider an amalgam of the second, third and 
fourth types of randomness composed of a "frozen" part 
(due to the mean positions of the network constituents) 
and a "molten" part (due to the thermal fluctuations of 
the network constituents). 



B. Overview 

In the present Paper we investigate the static structure 
of nematic alignment in IGNEs that can be probed in 
the high-temperature regime. This investigation is partly 
motivated by the experiments on polydomain structure 
reported in Ref. [9] . We focus our considerations here on 
systems in the high-temperature regime and which are 
not subject to externally applied deformations. 

Two types of nematic fluctuations are present in 
IGNEs: (i) those that are frozen in, however imperfectly, 
by the network during the process of cross-linking; and 
(ii) thermally driven departures away from the mean lo- 
cal alignment pattern that is frozen in. In order to char- 
acterize such fluctuations, we make use of (i) the correla- 
tor of the local nematic order that was frozen in during 
the cross-linking process, this correlator being appropri- 
ately averaged over realizations of the quenched disorder 
(and termed the glassy correlator); and (ii) the corre- 
lator of the thermally driven departures of the nematic 
order from the mean frozen-in local alignment pattern, 
appropriately averaged over realizations of the quenched 
disorder (and termed the thermal correlator). 

We have already mentioned that the physical proper- 
ties of IGNEs depend on the conditions under which they 
were prepared. To reflect this fact, we make a careful dis- 
tinction between two thermodynamic ensembles in our 
theoretical approach: the first, in which the IGNE was 
prepared, which we term the preparation ensemble; and 
the second, in which the system is measured, which we 
call the measurement ensemble. This structure enables 
our theory to capture the ability of IGNEs to "remember" 
at least to some degree, the local nematic alignment pat- 
tern at the moment of preparation, and in addition en- 
ables the determination of the dependence of the strength 
of this memorization on (i) the temperature at which the 
system was cross-linked, and (ii) the average number of 
cross-links per polymer. 

We note in passing that having both preparation and 
measurement ensembles places us in a family of dis- 
ordered systems for which spontaneous replica symme- 
try breaking is expected to be irrelevant, in contrast 
with settings that do not feature a preparation ensem- 
ble; cf. Ref. [20]. (Technically, this expectation shows up 
in the need to investigate coupled replicas of the physical 
system in the neighborhood of not zero but one replica.) 
Thus, we do not expect our approach to yield glassy phe- 
nomena such as hysteresis in the stress-strain behavior of 



IGNEs — which is not unreasonable, given the absence, to 
date, of experimental observations of such phenomena. 



The fact that the elastomer network is thermally fluc- 
tuating means that any prospective theory of IGNEs 
should feature a typical localization length-scale, below 
which the polymer constituents of the network are effec- 
tively delocalized. Our approach features such a length- 
scale, which leads to the possibility that nematic correla- 
tions undergo a novel, oscillatory form of decay with dis- 
tance in a certain regime. Prior theoretical approaches 
built on conventional random-field models do not fea- 
ture the thermal fluctuations of the elastomer network 
(see, e.g., Refs. [12, 13, 19-21]), and thus do not cap- 
ture this intriguing phenomenon. A phcnomcnological 
random-field-type model that does take the thermal fluc- 
tuations of the elastomer network into account was pre- 
sented in Ref. [22]. In the present Paper we derive the 
Landau-type free energy for IGNEs that was presented in 
Ref. [22], doing so via a microscopic model that involves 
dimers that are randomly and permanently connected by 
Hookean springs. The present Paper thus provides a mi- 
croscopic justification for the ideas and results presented 
in Ref. [22]. 



The outline of the Paper is as follows. In Sec. II we 
present a microscopic dimer-spring model of IGNEs. In 
Sec. Ill we apply the replica technique and implement 
the Hubbard-Stratonovich scheme to decouple the in- 
teracting microscopic degrees of freedom. In Sec. IV 
we derive a Landau- Wilson type of free energy for the 
IGNE, which involves an order parameter field Q for 
the isotropic-to-nematic phase transition as well as an 
order parameter field f2 for the vulcanization/gelation 
transition. In Sec. IV A we make an expansion of the 
Landau- Wilson free energy for small Q and J7, which 
is appropriate for exploring the physics in the vicinity of 
the transition to the random solid state. In Sec. V we de- 
termine the stationary states of the Landau- Wilson free 
energy. In Sec. VI we derive an effective replica Hamil- 
tonian describing local nematic order in IGNEs by set- 
ting 17 to its stationary value but retaining fluctuations 
of Q to quadratic order. In Sec. VII we then compare 
this effective Hamiltonian with that arising from the phc- 
nomcnological Landau free energy considered in Ref. [22] 
and show that they are equivalent. In Sec. VIII we use 
the effective Hamiltonian to derive the glassy and ther- 
mal correlators. In Sec. IX we describe two alternative 
microscopic models of the IGNE, viz. (i) a worm-like 
chain model of side-chain nematic polymer networks; and 
(ii) a jointed chain model of main-chain nematic polymer 
networks. As we shall see, at length-scales larger than 
the size of a nematogen, both of these models give rise 
to liquid-crystalline characteristics similar to those re- 
sulting from the dimer-and-springs model, reflecting the 
fact that the three models inhabit a common universality 
class. In Sec. X we make some concluding remarks. 



II. INGREDIENTS OF THE MODEL 

We model an IGNE microscopically as a system of N 
dimers in D spatial dimensions that are randomly, per- 
manently linked via springs (see Fig. 1 and Ref. [23]). 
We envision the springs as mimicking the flexible con- 
stituents of liquid crystalline polymers whilst also serv- 
ing as cross-links; the dimers mimic the stiff constituents 
of liquid crystalline polymers. Each dimer (labeled by j, 
where j = 1, . . . , N) consists of two particles at position 
vectors Cj^i and Cj^_i separated by a fixed distance i. 
The orientation of the j-th dimer is specified by the unit 
vector 

n, = ^^Z^l^f-L. (2.1) 

The dimers interact via three types of forces. First, there 
is an orientational interaction between dimers that pro- 
motes parallel or antiparallel alignment. We model this 
interaction via a potential of the Maier-Saupe type, viz., 

N 



H„ 



V 



^ J(c,-c,)(n,.n,)2, (2.2) 

where Cj = (cj.i + Cj^_i)/2 is the position of the j-th 
dimer's center of mass. We assume that the aligning in- 
teraction is short-ranged, and model the interaction po- 
tential J{c) by the form {Jo/{2na'^)°/^) expi-c^ /2a'^). 
In Fourier space, the potential is given by Jp — 
Jo exp(— p^a^/2), where a specifies the range of the inter- 
action between dimers and Jq characterizes its strength. 
In addition to the orientational interaction, the dimers 
experience a positional excluded-volume interaction be- 
tween particles belonging to any pair of dimers, which we 
model via an Edwards- type pseudo-potential [24, 25] : 
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(2.3) 



where A is the strength of the excluded-volume interac- 
tion [26]. The presence of sufficiently strong excluded- 
volume forces stabilizes the system against collapse to a 
globule, even when well cross-linked. Third, any given 
two dimers that are connected by a spring are taken to 
interact additionally via a harmonic potential associated 
with the spring, which we take to have zero rest-length 
and native mean-square end separation b^ , characteristic 
of Gaussian molecular chains. These springs and, specifi- 
cally, the architectural information indicating which pairs 
of rod ends are connected to one another by springs con- 
stitutes the quenched randomness x oi any given real- 
ization of the system. This information takes the form 
X = {ie,Se,je,te}eLi, wherc M is the total number of 
springs and spring e connects end Se of rod i^ to end ig 
of rod je ■ These springs result in the following Hookean 
term in the Hamiltonian: 
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FIG. 1: A model of dimers (indicated by dumbbells) cross- 
linked via harmonic springs. The dimers have length £ and 
the springs have r.m.s. length b. 



where T denotes the temperature and we have adopted 
units in which Boltzmann's constant is unity. 

The total Hamiltonian for the dimcrs-and-springs 
model is then given by 
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(2.5) 



For a given realization x of th^ quenched disorder, the 
corresponding partition function takes the form 



Zx = 



/ ]Jdci_idci_i(5(|cj_i -Cj-i| -£) I 



-H^/T 



and the free energy F^ is given by —TlnZ^. As is well 
known [27], it is appropriate to average the free energy 
of the system over realizations of the quenched disorder. 
Denoting this average by the square brackets [•••], one 
has 



[^x]=E^x^x--rE^xln^, 



(2.7) 



III. REPLICAS AND COLLECTIVE FIELDS 
A. Statistics of quenched disorder 

What statistical distribution should one use to com- 
pute the average of the free energy over the quenched 
disorder? In common with other elastomers such as 
isotropic rubbery systems, the quenched-disorder average 
for IGNEs can be performed via a variant of the Deam- 
Edwards distribution P{x) [28]. Such distributions re- 
flect situations in which systems undergo instantaneous 
cross-linking: one begins with a melt or solution at equi- 
librium and — so rapidly that hardly any relaxation has 
time to occur — one introduces permanent bonds between 
some random fraction of the pairs of dimers that hap- 
pen, at the instant of cross-linking, to be nearby one 
another. To construct the associated P{x)i one should 



respect the causal order of the cross-Hnking process (see, 
e.g. Ref. [30]). FormaUy, this amounts to the specifica- 
tion 



Pix)^JVc"P{x\{cl})Pi{cl}), (3.1) 

where the multiple integral J Dc^ is defined via 

j Vc'' ^f[ j dcl_,dcl,, (3.2) 



the distribution P{{c^ g}) describes the statistics of the 
liquid system at the instant of cross-linking, so that 
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and ^(xllc",,}) is the conditional probability that cross- quircmcnt that ^(x|{c"g}) be properly normalized over 
links arc formed between pairs of dimer-ends at positions y j.e. V] ^(xljc" }) = 1' it is given by 
{(c^__ g^, c°^ tJ}^=l^ given that the constituents of the liq- 
uid at the instant of cross-linking are at {c^ g]- It is given 
by [30] 
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(3.5) 
(3.4) The probability P{x\{(^^ s\) features a dimensionless pa- 
rameter r] that controls the likelihood that cross-links are 
actually formed. In App. A we show that it is related to 
Here, V is the dimensionless volume of the system, i.e., the average number of cross- linking springs per dimer via 
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y/(27r62)^/2. The term ffnorm(c°) arises from the re- the formula [M]/N « 2ff 



B. Disorder averaging the free energy 

Now that we have constructed a suitable disorder distribution P^, we use it to perform the disorder average of 
the logarithm of the partition function [InZ^], doing this indirectly, using the replica technique (see, e.g., Ref. [31]). 
Thus, we represent the logarithm in Eq. (2.7) as a limit to obtain 
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where we have interchanged the order of taking the replica limit and performing the disorder average on going from 
Eq. (2.7) to Eq. (3.6). We then insert the Deam-Edwards distribution (3.1) to obtain 
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Here, a — 0, 1, . . . , n labels the replicas, and Zf? is dc- Note that in Eq. (3.7) the term -ffnormdc"^}) limits per- 
fined via mutation invariance to replicas a = 1, . . . , n, with replica 

N 

xTT^dc" -cO 1-^). (3.8) 



a — being excluded. This is consistent with our phys- 
ical expectation that the preparation and measurement 
ensembles play distinct roles. 

Next, we complete the computation of the disorder 
average in Eq. (3.7), which yields 
replica partition function ^i+n, i.e. 



Z"] in terms of the 
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in which J dc = Y[d=i J dcd, where the effective replica 
Haniiltonian iJrcp is given by 
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We note that in iJrep the replicas are coupled as a result of the disorder averaging. See App. B for details of the 
derivation of Eqs. (3.10,3.11). 



C. Collective fields and their physical meaning 

We note that in ifrcp, Eq- (3.11), the N interact- 
ing (replicated) dimers are coupled with one another, 
which means that the trace over microscopic variables 
Ci^s in ^i+n [Eq. (3.10)] cannot be straightforwardly car- 
ried out. We decouple these dimers by defining micro- 
scopic collective fields cof, and qtp, and making a Hubbard- 
Stratonovich transformation involving fiuctuating auxil- 
iary fields conjugate to the collective fields. Thus, we 
arrive at a description in terms of N uncoupled copies 
of a single replicated dimer, for which the trace over Ci,s 
can be readily performed, order by order in an expansion 
in the auxiliary fields. The collective fields are given by 
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Here, we have used the symbol k to denote the n + 1- 
fold replicated wave- vector (fe°, k^, . . . , fc"), and have re- 



stricted the value of k in oj^ to the "higher-replica sector" 
(HRS), viz., the set of replicated vectors each having at 
least two non-zero vector entries. Details of the Hubbard- 
Stratonovich transformation are given in App. B. The re- 
sult is that the replica partition function Zi+„ becomes 
a functional integral over the auxiliary fields Q^ and Qp, 
whose expectation values are related to the expectation 
values of luj^ and Qp via 



(Qp)/ 






(3.13a) 
(3.13b) 



as we demonstrate in App. C. Note that (•••)/ indicate 
an expectation value taken with respect to the Landau- 
Wilson free energy /i+„ discussed in Sec. IV and defined 
in Eq. (4.2), and {■■■)x denotes an expectation value 
taken with respect to the Hamiltonian H^, Eq. (2.5). 



To gain some intuition for the physical significance of [(f^)/] and [(Q")j], we transform Eqs. (3.13) to replicated 



real space: 
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Thus, except for a trivial constant, we can interpret 
{fl(f))f as the joint probability that a given dimer end 
is found at position r^ at the instant of cross-linking, 
and that the same dimer end would be found at n sub- 
sequent widely-separated time instants at the positions 
r^,...,r" [23]. {il{f))f vanishes if all dimers are delo- 
calized, and has a nonzero value if a fraction of them are 
localized. Thus, (0(f))/ serves as the order parameter 
that detects the phase transition from the liquid state 
to the random solid state. Similarly, (Q°(r))/ is the ne- 
matic order parameter for the preparation state, whilst 
(Q"(r))/ (for a = 1, • . . , N) is the nematic order param- 
eter for the measurement state. 

We note that when defining Q we have chosen to ex- 
clude the components associated with the one-replica sec- 
tor (denoted "IRS"), i.e., the set of replicated wave- 
vectors that each have only one non-zero vector entry 
[i.e., k = {0, . . . , A;", . . . , 0)]. This is because this sector 
of the field corresponds to fluctuations in the macroscopic 
density of dimers, and is strongly stabilized as a result of 
the excluded volume interactions, regardless of the extent 
of the cross-linking. The physical content of the decom- 
position of fields into higher and lower replica- sectors 
is that if condensation occurs in the higher sector only, 
this implies the random localization of particles. On the 
other hand, if condensation occurs in the IRS as well, 
this indicates the formation of a state having some kind 
of spatially modulated density structure. We are focus- 
ing on highly incompressible systems, for which density 
fluctuations are negligibly small. Incompressibility is en- 
forced by taking A to have a large value, so that fluctu- 
ations of IRS counterpart to Q are strongly suppressed. 
We also define the zero-replica sector (denoted "ORS") 
to be the set whose only member is the replicated wave- 
vector that has zero for every entry. The lower-replica 
sector (denoted "LRS" ) would then refer to the union of 
the one-replica sector and the zero-replica sector. 



D. Intermezzo on replicas 

Compared with the more familiar example of spin 
glasses (see, e.g., Ref. [31]), the effective replica Hamilto- 
nian (3.11) in our theory contains an extra replica. What 
meaning can we ascribe to this "zeroth" replica as well 
as to the other n replicas? Physically, the zeroth replica, 
which originates in the Deam-Edwards distribution, cor- 



responds to the state of the system at the instant of 
preparation, whilst the other n replicas correspond to the 
state of the system when it is measured. The dependence 
of the measured properties on the state of the system at 
preparation corresponds, operationally, to the coupling 
in the effective replica Hamiltonian (3.11) between the 
freedoms belonging to the zeroth replica and those be- 
longing to the other n replicas. On the other hand, such 
a coupling does not imply that the measured properties 
influence the preparation state, as one can show via a 
careful consideration of the n — > limit [30] . 

Because the replica framework accounts for both the 
preparation and measurement ensembles vis-a-vis the ze- 
roth and other n replicas, it is well adapted to investi- 
gations of nematic elastomers, whose measured proper- 
ties are known to depend on their preparation histories; 
see, e.g., Ref. [9]. Such a dependence was already un- 
derstood by Deam and Edwards and by de Gennes more 
than thirty years ago (see, e.g., Refs. [28, 33]). 



IV. LANDAU- WILSON FREE ENERGY 

As we demonstrate in App. B, the local incompress- 
ibility of IGNEs allows us to express the effective replica 
theory in terms of the auxiliary fields Q and Q" as 

/n 
Pr!^[]P0^exp(-7V/i+„(r!,Q")), (4.1) 

where the Landau- Wilson free energy per dimer /i+n is 
given by 

~2 1 " T 
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We have introduced the notation 



Gi{n) 



^"2 



2V 



-E^fe^E 



-ik-Cs 



(4.3) 



G2{Q) ^ EE^Q2,..(p)e-^^-^° 
p 

(nln2, - ldd,d?j , (4.4) 



a—Q p 

X 



and have specialized to three spatial dimensions (i.e., Moreover, J^k denotes the sum over replicated wave- 
D = 3). The curly braces {SS'} in Eq. (4.2) denote vectors that are restricted to the higher-replica sector, 
the trace of the product of the tensors S and S', i.e., 

E?i,d2=i Sd^d^S'd^d,- The symbols 77^ A^ and (• • Oia+ri 
are defined via 

if = Aff, (4.5a) 

Afe = exp(-fe2fcV2), (4.5b) 

A. Expanding the Landau- Wilson free energy 

To develop the expansion of the Landau- Wilson free energy we expand the log-trace term in Eq. (4.2) in powers of 
Q and O to obtain (see App. D for details): 

/i+„(f7, Q) « fnm + /q(Q) + ,fc{n, Q). (4.6) 

The first term on the right hand side, fii{n), describes the vulcanization/random solidification transition for an 
isotropic elastomer [29, 34, 37] and is given by 

k 

It exhibits a linear instability at the critical value 77^ = 1, [36] which reflects the destabilization of the liquid state with 
respect to a gel/random solid state when the average number of cross-links per dimcr is increased beyond a certain 
critical value. As IGNEs are random solids, they must have undergone a vulcanization transition (controlled by the 
density of cross-links) from the liquid state to the random solid state. Thus, we keep the J7-only terms to cubic order; 
this is in line with the approach adopted in Ref. [34]. The second term, /g, is given by 

n 

/q(Q) = EE^^'^(^"'P){QpQ-p}' (4-8) 
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These terms describe the free-energy cost of inducing nematic alignment from the unaligned state. As the focus of the 
present Paper is on the liquid-crystalline properties of IGNEs at high preparation and measurement temperatures, 
we retain the Q-only terms to quadratic order. This collection of terms is an (1 + n)-fold replicated version of 
the (quadratic part of the) Landau-de Gennes free energy for a liquid of nematogens at high temperatures [35] . In 
Eq. (4.8), r° is the preparation temperature (i.e., the temperature at which the system was cross-linked), and T" 
(for a = 1, . . . ,ri) is the measurement temperature (i.e., the temperature at which the system is measured, long after 
cross-linking). As replicas a = 1, . . . ,n are actually copies of a single system measured at one temperature, we have 
T" = T for a = 1, . . . , n. The kernel 7W(T",p) is given by 

which, expanded in powers of wave-vectors, yields [32] 

Al(T",p)«yt"r + £V, (4.10) 

where t° is the reduced preparation temperature (T'^ — T*)/T* (where T* has the value 4Jo/15), and i" (for a — 
1, . . . , n) is reduced measurement temperature, (T — T*)/T* . For simplicity, we write t" = t for a = 1, . . . , n. The 
parameter A" (= J^T* /T") characterizes the aligning tendencies of nematogens, and £" [= (2T* — T°')a^ Jq/2T"] is 
the generalized nematic stiffness in the (Landau-de Gennes equivalent of the) one- Frank-constant approximation [35] . 
Hence, the nematic free energy /g becomes 

1 " 1 

/q(Q) « ^ ^ M(T", p){Q^ Q%} + E E ^(-^^ + ^P'){Qp Q-p}' (4.11) 

p a—1 p 



where A and C are, respectively, A" and C" ioi a — 1, . . . ,n. The third term of the right hand side of Eq. (4.6), fc, 
describes the couphng between the nematic freedoms and the elastomer and is given by 

n 

a^j3 p,q 

n 

n 

+ E EEr^frf.<i3rf.(P''?'^)^fe"-^-pe--..^^Q?.<i.(P)<d.('7) 

n 

+ E E E \^d,d2d:,dAP^ ^) Vq.O + ^Xd2d-sdSP^ 9' ^)) 

X f^i ^-fe-(p+,).- QldM Qldjci), (4.12) 

where the wave- vector-dependent coefficients $, ^, 8, T and F are defined as 

~2 T T 

*dfd2d3d4(P'9) = -A_pgc_qgg ;^6ooya ^p PdiPd^qd.qd,] (4.13a) 

(4.13b) 

T?,d.d3d.b,fc) ^ -A|^(^^) <5rf,rf3<5rf,d,; (4.13c) 

f^ifi J J 

'^ did2d3d4\P>'i''^) — ^^QrparpP ^k^-k-pe^-qil^ 

X (fc." fc." fcl^l + J (P.^fc.. + fc." p..) (^.afcl + fcl?.. 

+ oP<iiPd2Qd3qdi ] ; (4.13d) 



■^.1.2.3.4^^' 9''^^ — 3QCyap^fe^-fe-(p+q)e° 



^2|p + q|2 t\p + q\'^ P\p + q + Ak"^^ 



56 8064 56 

t\p + q + Ak°'\^ 



8064 

1 ^2|p^q|2 



(5.1 .3(5^2^4 

^^(P + '?)<il(P + 9)d3 '5.2.4 
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1 £2|p^q^4^a|2 

14 ^ 1008 

x^2(p + g + 4fc")d, (p + q + 4fc")rf3(5rf,d4 

'^ 2016 "-^ + 9).l (P + <?).2 [P + 'Z)d3 (P + 9).4 

+ 2^6^^ + '^ + ^^"^"^^^ + '^ + '^^"^''^ 

X (p + g + 4fc")<i3(p + q + 4fc")dJ . (4.13e) 

Here, we have retained the lowest-order terms that couple fi and Q, i.e., those proportional to QQJ7 and Qf7f7. 
We have, moreover, kept only the lowest-order wave-vector dependencies of their coefficients. We have, in addition. 



retained the term proportional to Q Q fi fi because we have found that it is the leading one responsible for the 
qualitatively new possibility of thermal and glassy correlators that not only decay but also oscillate with distance. 
The physical significance of fc is that it encodes into the theory the fact that the ncmatic degrees of freedom 
inhabit an environment that is at the microscopic level anisotropic, inhomogcncous and thermally fluctuating, but 
that no remnants of this anisotropy or inhomogeneity survive to the macroscopic level, fc thus also describes the 
mutual effects of random localization and ncmatic alignment. In particular, it enables us to study the impact of the 
memorization of the ncmatic fluctuations in the preparation state on the ncmatic fluctuations in the measurement 
state [22]. 



V. SADDLE-POINT APPROXIMATION 

Next, we consider the saddle-point equations for Q and 
n, which follow from the stationarity of Eq. (4.6): 

Sfl+n 




0; 



1 r rp Sfl+n 



(5.1a) 
= 0. (5.1b) 



As we are concerned with macroscopically isotropic 
states, we require that the saddle-point value of Q van- 
ishes: Q = 0. As a result, the saddle-point equa- 
tion (5.1a) for f2 is the same as it would be for ordinary 
(i.e., non-nematogenic) elastomers, and reads 



in which we have defined the length-scale £'j. via £'^ = b^ + 
(£^/6). Next, by defining the quantities e = 3(?7^ — l)/ri* 
and ir = VStj./rj'^, the saddle-point equation takes the 
form analyzed in Ref [37]: 

^ + bl\k\')n^-^Y.^yn;_y = o, (5.3) 



This equation is automatically satisfied, which we can see 
as follows. The product ^fc^_fe_„ea involves a product 
of two Kronecker deltas (5j^"_ fc^.o'^-p-X]"- fc°,0 7 which 
reflects the macroscopic translational invariance of the 
random solid, and implies that p — 0. Next, we note 
that the restricted wave-vector sum ^^, can be replaced 
by the full wave- vector sum '^f., as the IRS and ORS con- 
tributions to the restricted sum vanish. Lastly, 
a scalar, implying that 



%^- 



T.^k^-kKki-l\nH,,a. 



= 0. 



(5.8) 



For simplicity, instead of working with the full distribu- 
tion P{t), we replace it by one that is sharply peaked at 
1/Ci- This approximation is valid as long as we are con- 
cerned with the broad implications for the nematogens of 
the very presence of localized network constituents. We 
expect that the spread of localization lengths would at 
most result in the quantitative but not qualitative mod- 
iflcation of our results. If the localization is sharp, the 
Ansatz for the order parameter becomes 



^k 



G 



dz 
— ( 

V 



'Ec'^^-^-jIfel'?? 



(5.9) 



and is therefore solved by mean of the the Ansatz 



% = G 



dz 



dTP(T)e'2:„fc" — I&IV2T 



(5.4) 



(with fc G HRS), provided that the gel fraction obeys 

G ^ 2e/3, (5.5) 

and P{t), which has the meaning of the distribution of 
inverse square localization lengths r, obeys 

_ (5.6) 
The typical value ^^ of the localization length is e ^^"^ir, 
which diverges as the vulcanization transition is ap- 
proached from the solid side. The saddle-point equation 
for Q", viz., Eq. (5.1b), yields 



E%^- 



k—pe^ 



[{pd.Pd, + (fc" +p/2)di(fc" +p/2)d,) 

(5.7) 



--(|p|2 + |fe"-Hp/2|)<5,,rf,j=0. 



VI. EFFECTIVE THEORY OF STRUCTURE 
AND CORRELATIONS IN IGNES 

In App. E, we derive an effective replica Hamiltonian 
for the liquid crystalline behavior of IGNEs, using an ap- 
proximation in which fl is set to its saddle-point value Cl 
in the Landau- Wilson free energy (4.6), but allowing Q" 
to undergo fluctuations. By making this approximation 
we are neglecting the impact of fluctuations in 51 on the 
nematic freedoms. We adopt this level of description be- 
cause it is the least complicated one that is capable of re- 
vealing the impact of the localized network on the liquid 
crystallinity characteristic of IGNEs. As it remains con- 
stant, we do not need to consider the fa contribution to 
the free energy. Furthermore, at the present level of ap- 
proximation, the term in fc proportional to fifiQ" van- 
ishes, as we show in App. E, as a result of the macroscopic 
translational and translational invariance of the random 
solid state. Moreover, as we show in the same appendix, 
at wavelengths long compared with ^l the dominant con- 
tribution to fc comes from a term in fi fj Q" Q^. Thus, 
in forming an the effective replica Hamiltonian for the 
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high-temperature hquid crystaUinity of IGNEs, which we 
denote by Hi+n[{Q°'{-)}2=o]' ^^ need only consider the 
/q term together with a term in fj fj Q" Q^ term: 



n 

a—l p 
1 "" 



a — l p 



1 " 

■^ E E^p{QpQ-p}- (6-1) 



,/3 = l p 



VII. PHENOMENOLOGICAL CONTENT OF 
THE MICROSCOPIC REPLICA THEORY 



We now show that the effective rephca Hamiho- 
nian (6.3) that we have derived on the basis of the mi- 
croscopic dimers-and-springs model can be interpreted 
as having arisen from a more phenomenological, contin- 
uum description of a liquid crystalline systems subject 
to a novel pair of interrelated random fields. This phe- 
nomenological description and its implications were ex- 
plore in Ref. [22]. Our purpose here is not to revisit 
these issues in detail but rather to reveal the microscopic 
underpinnings of the phenomenological theory. In that 
latter theory, the free energy of an IGNE having a given 
realization x' of the quenched disorder is given by 



Let us draw attention to the kernel Hp, which is defined 
via 



Hp =^ Hqc 



2P 4l 



(6.2) 



Here, Hq, which we call the disorder strength, char- 
acterizes the strength with which the network in- 
fluences the liquid crystaUinity, and has the value 
G^Jifi^£^/{900T^l). The kernel Hp is the manifesta- 
tion of the presence of the thermally fluctuating random 
elastomeric network and, in particular, encodes the cen- 
tral physical characteristic of the network, viz., that at 
long length-scales the localization appears perfect but at 
length-scales shorter than ^^ thermal fluctuations render 
the network effectively "molten'.' Note that the charac- 
teristic length-scale beyond which H{r) is suppressed is 

a- 

To determine the effect that the preparation history 
has on the equilibrium liquid crystalline properties post- 
cross- linking, we integrate out the zeroth- replica element, 
Q°, thus obtaining the effective Hamiltonian 



(At + Cp^ 



a—l p 



\Hr.\' 



\ 2T 



H, 



2T0{AH0+C0p2^ 
\HJ^ 



{QpQ-p} 



5i ^ \2T + 2rO(yt/+£V) j ^^P^-P^-^^-^^ 

A noteworthy feature of iJ^ff is that the replica-diagonal 
contribution to the quadratic term in Q" is structurally 
distinct from the replica off-diagonal contribution. This 
structural feature, as well as the short length-scale liq- 
uidity encoded in Hp, enables us to capture a richer 
range of physical behavior (such as oscillatory-decaying 
nematic correlations) than can be predicted via conven- 
tional random-field approaches, for which the replica- 
diagonal and off-diagonal terms contain identical coef- 
ficients. 



^x'[Q] = \Y,({At + Cp' + Hp){QpQ^p] 



-2{{Y p + {T/T^)HpQl)Q_p]), 



(7.1) 

where x' = {Q°jY}, and Y and Q" are independent, 
Gaussian-distributed random fields, with zero means and 
non-zero variances, the latter being given by 



3p+p',0 



[jYpYp'IJ = T HpSpj^piQ . 



(7.2a) 
(7.2b) 



Here, Q*^ describes the impact of the configuration of the 
nematogens that is present at the instant of cross-linking 
on the post-cross-linking nematic alignment pattern, and 
Y accounts for the impact of the local anisotropic en- 
vironment created by the localized network constituents 
post-cross-linking on this nematic alignment pattern. For 
a given realization of x', the partition function is given 

by 



Z^, ^ J VQeM-H^'[Q]/T), 



(7.3) 



and by suitably averaging over x' using the replica tech- 
nique we obtain 



[^x'] 



ex 



/n / n ^ 

n2?Q"exp -E^x'[Q"]/r 

-J []PQ"exp(-i7'[Q"]), 



(7.4) 
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where the effective repHca Hamiltonian i?phcn[Q"] is 
given by 

-ffphen[Q"] 

1 " 



2T 



Q — 1 p 



1-2 E E[iYpY_p}]{Q^Q^p} 



2T2 



a, 13=1 p 



^ E Ei^pn{Q°Q°}]{QpQ-P 



}■ 

(7.5) 



Applying the disorder statistics specified in Ref. (7.2), 
we arrive at the result that i?phcn[Q"] = -ffcff[Q"], i-e., 
the phcnomcnological continuum description originally 
reported in Ref. [22] contains the same physics as the 
microscopic dimers-and-springs model. 



VIII. STRUCTURE AND CORRELATIONS IN 
IGNES 

A. Diagnostic quantities 

To describe the essential features of the pattern ne- 
matic ordering characteristic of IGNEs in the high- 
temperature regime we focus on the following pair of 
thermally- and disorder-averaged correlators: 
(i) the average, taken over realizations of the cross- 
linking, of the product of the local nematic order at two 
points, i.e.. 



C^(r,r')^[{(q(r))(g(r')})] 



(8.1a) 



which we term the "glassy correlator" ; and 

(ii) the disorder average of the familiar correlator of the 

fluctuations in the nematic order, i.e., 

C^(r,r') ^ [({(g(r) - (g(r))) (q(r') - (q(r')))})], 

(8.1b) 
which we term the "thermal correlator." 

The correlator C^ characterizes the strength of the 
thermal fluctuations of the nematic alignment away 
from the local mean value as well as the spatial range 
over which these fluctuations are correlated. Inter alia, 
through its range, C^ is capable of signaling the occur- 
rence of a continuous phase transition. The correlator 
C^ is a diagnostic of particular value for nematic elas- 
tomers, as it detects the occurrence of randomly frozen 
(i.e., time-persistent) local nematic order. For the case 
where r and r' are co-located, it is the nematic analog 
of the Edwards-Anderson order parameter, introduced 
long ago for spin glasses [38], in the sense that it mea- 
sures the magnitude of local nematic ordering, regard- 
less of the orientation of that ordering. Moreover, how 
C varies with the separation of r and r' determines the 




FIG. 2: Schematic depictions of snapshots of nematogen lo- 
cations and orientations at a particular instant (blue, full), 
and at a much earlier instant (gray, shaded), (a) A conven- 
tional liquid crystal in the isotropic state. Such systems do 
not memorize the local pattern of nematogen alignment in- 
definitely. There is no correlation between the orientations 
of blue and shaded nematogens that are depicted near one 
another. Nor is there any preference for blue and shaded ne- 
matogens that are depicted near one another to be the same 
nematogen. (b) A liquid cystalline elastomer in the macro- 
scopically isotropic state. Such systems do memorize the local 
pattern of nematogen alignment indefinitely. The orientations 
of blue and gray nematogens depicted near one another are 
likely to be correlated. For systems in which the nemato- 
gens are chemically bonded to an elastomer network, blue 
and shaded nematogens depicted near one another are likely 
to be the same nematogen. However, for systems in which the 
nematogens are not chemically bonded to a network there is 
no preference for blue and shaded nematogens depicted near 
one another to be the same nematogen. 



spatial extent of regions that share a roughly common 
nematic alignment. Two mechanisms are responsible for 
the existence of these aligned regions. First, the forma- 
tion of a random network causes a local breaking of ro- 
tational invariance, which has the effect of creating ran- 
domly anisotropic environments that tend to align the 
nematogens locally. Second, although the equilibrium 
state of the system at the instant prior to cross-linking 
is, on average, isotropic, a "snapshot" of its microscopic 
configuration at that instant would reveal local nematic 
order of the type that we normally call thermal fluctu- 
ations. The cross-linking process can trap these fluctu- 
ations in, either partially or fully, the extent depending 
on the strength of the cross-linking and the temperature 
at the moment of cross-linking. 

The correlators C^ and C'^ have been computed in 
Ref. [22] via Eq. (7.1). In the present section, we re- 
derive the results via the replica approach to the micro- 
scopic dimers-and-springs model that we have developed 
in the foregoing sections. To proceed, we make use of 
the effective replica Hamiltonian (6.3) and the following 
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identities (which we prove in App. C): 

[{(g(r))(9(r'))}]=nm((Q"(r)Q^(r'))) 

n— J-0 

{a^/3), (8.2a) 
[{{q{r)q{r'))}] = lim ((Q"(r)Q"(r'))), (8.2b) 

n— fO 

where ((. . .)) denotes an average performed with respect 
to the effective Hamiltonian H^g, Eq. (6.3). To compute 
the quantities on the right hand side of Eqs. (8.2), we 
invoke the quadratic form of Eq. (6.3), invert the cor- 
responding kernel, and use the rephca diagonal and off- 
diagonal parts to obtain 



a =5T 



5T 



1 



At + Cp^ + iJj 



2^0 \-^ ^ 



P 

£V)" 



'\Hp\' 



H„ 



{At+Cp^ + Hp 



(8.3a) 
(8.3b) 



Note the presence of the scale-dependent kernel Hp in the 
denominators of the correlators, which plays an essen- 
tial role in determining their behavior. This may appear 
surprising when we compare with the result one would 
obtain via a conventional random-field approach (for de- 
tails, see App. F). At length-scales larger than ^l, the 
presence of Hp in the denominator of Cp leads to a down- 
ward renormalization of the bare critical temperature T* 
by an amount proportional to the disorder strength (and 
hence grows with the cross-linking density). This indi- 
cates that the nematogens are more strongly inhibited 
from aligning with one another if the density of cross- 
links is higher. Interestingly enough, a similar result has 
been obtained using the molecular level neo-classical elas- 
ticity theory of nematic elastomers [17]. This disordering 
effect of the random polymer network (on the nematic 
alignment) is not contained in the conventional random- 
field type of models [40]. In addition, a larger value of 
Hq leads to a larger amplitude of Cp , indicating that 
a stronger localization of the network results in a more 
strongly trapped-in nematic pattern. Also note that the 
correlators Cj and Cp meet the physically sensible re- 
quirement that they revert to the forms appropriate to 
a nematic liquid in the absence of a network at length- 
scales very short compared with ^i, for which Hp be- 
comes very small, indicative of the molten character of 
the network at such scales. 



B. Oscillatory-decaying correlations 

The scale-dependent kernel Hp in the denominator of 
each of the correlators in Eq. (8.3) also gives rise to 
the possibility that the correlators undergo both oscil- 
lation and decay with distance. As discussed in Ref. [22], 
C'^{r) undergoes oscillatory decay with distance for Hq > 
i/'^^''(= 2£/^|), and the length-scale of the oscillation is 
given by £,t,o 
dent of T. 



^l/V^^^Jh^JW^, which is indepen- 
On the other hand, C'-^{r) undergoes oscil- 




C^(r) 



(b) 



■r/& 



glassy 



FIG. 3: Real-space decay behavior of (a) the 
relator (rescaled) C^(r) = {12tv^ C/T)C'^{r), for i° > 
THo/T°A°; t = O.lC/Ail, at (i) Ho/H^"'' = 0.5 (weak 
disorder; red, dashed) and (ii) Ho/H^'^' — 40 (strong dis- 
order; blue, solid), (b) the thermal correlator (rescaled) 
C"^ {r) = (27r^£/5r) C"^(r), for the same parameters. Ongoing 
from weak to strong disorder, both correlators cross over from 
simple exponential decay to oscillatory decay of wavelength 
of order ^l ■ 



H 

6 



Thermal and 
glassy oscillations 



Thermal oscillations 



No oscillations 

2 



FIG. 4: Crossover diagram for the glassy and thermal cor- 
relators, indicating the three qualitatively distinct regimes of 
behavior for a system that is cross-linked at a very high tem- 
perature. Here, H = Ho/ H^ is a measure of the disorder 
strength, and t is the rescaled reduced temperature, with the 
value {A^\/ C)t. Above the blue solid line, both correlators 
oscillate and decay as a function of separation. Between the 
blue solid and red dashed lines, both correlators decay but 
only the thermal one also oscillates. Below the red dashed 
line, both correlators decay but neither oscillates. 



latory decay for sufficiently large values of Hq and suf- 
ficiently low values of T (see Fig. 3), with its oscillation 
scale £^G,o given implicitly by the following equation: 

1 + {iN/iG.of + ^{S,N/S.Lf - (i/oMt)e-«' /2«o.° = 0. 

The cross-over boundary between the oscillatory and 
non-oscillatory regimes for C"-^, which is the threshold at 
which the inverse oscillation wavelength I/^g.q, increases 
from zero to a non-zero value, is given by 



Ho=At + 2H 



(c) 



(8.4) 



The different regimes of behavior of C"^(r) and C^{r) are 
shown in Fig. 4. 



13 



IX. ALTERNATIVE MICROSCOPIC MODELS 
OF IGNES 

Instead of the dimers-and-springs model of IGNEs, one 
could have started with alternative microscopic models 
that may at first sight appear to be more faithful repre- 
sentation genuine IGNEs. For example, one could have 
started with either of the following two microscopic mod- 
els: (i) Model A: a system comprising worm-like chains — 
pairs of which are permanently, randomly bonded by 
point-like cross-links — as well as stiff rods that dangle 
fromn each chain at regular arc-length intervals; and 
(ii) Model B: a system of chains, each constructed from 
stiff rods that are connected in series and then perma- 
nently bonded by point-like cross- links between randomly 
chosen rod ends. Models A and B are, respectively, car- 
icatures of side- and main-chain nematic polymer net- 
works; starting from either Model A or Model B, we can 
apply the Hubbard-Stratonovich decoupling procedure 
described in Sees. II to IV, and thus derive a correspond- 
ing Landau- Wilson free energy that is structurally equiv- 
alent (i.e., having the same symmetries and types of cou- 
plings) to that derived for the dimers-and-springs model; 
cf. Eq. (4.2). The purpose of this section is to estabhsh 
this structural equivalence. We can then, in principle, 
compute the coefficients of the terms of the expansion of 
the Landau- Wilson free energy separately for Models A 
and B, apply a similar saddle-point analysis, and derive 
an effective replica Hamiltonian that would enable us to 
explore nematic fluctuations. However, the coefficients 
of the terms of such expansions are technically more dif- 
ficult to compute than those for the dimers-and-springs 
model. Moreover, even if one were to succeed in com- 
puting such coefficients, the expansion would still lead, 
at sufficiently large length-scales, to predictions that are 
identical, up to an overall length-scale, to those made 
on the basis of the dimers-and-springs model, Eq. (4.6). 
This is why we chose to work with the dimers-and-springs 
model in deriving an effective Hamiltonian for the liquid 
crystallinity of IGNEs. 

In what follows, we consider Models A and B sepa- 
rately. We define the microscopic Hamiltonian, introduce 
the corresponding collective fields, perforin the Hubbard- 
Stratonovich decoupling procedure and, lastly, carry out 
the log-trace expansion. We then show that the auxiliary 
fields and the terms in the Landau- Wilson free energy 
corresponding to either model have the same symmetries 
and structure as those obtained from the dimers-and- 



springs model, and thus, give rise to effective Hamiltoni- 
ans that are structurally equivalent to Hi+n] cf. Eq. (6.1). 
To streamline the presentation we display only the essen- 
tial equations, as the formal procedure employed in this 
section is identical to the one employed in Sees. II and 
III. 
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FIG. 5: Model of a side-chain nematic polymer consisting of 
nematogens (red online) and worm-like chain backbone (blue 
online). The chains are labeled by i (= 1, • ■ • , -'V), is of length 
L, and has S nematogens attached. The segment at arc-length 
T of chain % has position Ci(j)- The nematogens are attached 
to the backbone at equally spaced arc-length intervals. 



A. Model A: side-chain nematic polymer netw^ork 

We represent a side-chain nematic polymer as a worm- 
like chain (see, e.g., Ref. [39]) of length L with nemato- 
gens attached at equal intervals of arc-length LjS along 
the chain. Each nematogen is represented by a unit 



vector rii » , where i 



1, 



, N labels the chains, and 



s — 1, . . . , 5 labels the nematogens. We denote the arc- 
length measured from one end of the chain by r (with 
< r < L). A segment at arc-length r along chain i has 
a position vector Ci(T)\ the s-th nematogen resides at 
arc-length sL/ S^ and its base, which is attached to chain 
i, has a position vector Ci{sL/S); see Fig. 5. We require 
that nematogens and tangent vectors to the chains in- 
teract via orientational forces that favor parallel or anti- 
parallel alignment. 

The side-chain nematic polymer network is formed by 
randomly, permanently cross-linking M pairs of such 
chains via point-like cross-links. The pairs of cross-linked 
chains and the number M are fixed for a given realization 
of quenched disorder, but vary across such realizations of 
it. In the absence of cross-links the system is specified 
by the following Hamiltonian: 
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ij = l 
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^ ^ J(c,(sL/^)-c,(iL/5))K,.n,-,)' 



-i.j — 1 s,i— 1 

A^ 5 



A 1 ^ 1 "^ 



i,j — l s—1 



\ dr 



"J.s 



(9.1) 



The partition function corresponding to this system, sub- 
ject to M cross-hnking constraints, is given by 



Zsc « II fvc,f[ f dn, 

i=l'' 3 = 1-' 

M 



Hsc/T 



(9.2) 



Here, J T>c = YIo<t<l I '^^i''')' ^^'^ I ^^ denotes in- 
tegration over the unit sphere of nematogen orienta- 
tions. To proceed further, we can follow the strat- 
egy described in Sees. II and III, including making a 
Hubbard-Stratonovich transformation to introduce the 
auxiliary fields {^j, f^" Q"}, to arrive at the correspond- 
ing Landau- Wilson free energy /gc: 



/SO = |^E,l^^l^ + ^EE^?l^."l^ + ^EE^Tr|Q 
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(9.3) 



r 



Here, A" = A — ri'^{6"''^ — l)/(iVy"^^), reflecting the state. Note that the Boltzmann average 
fact that the excluded-volume parameter A is not renor- via 
malized by the effects of cross-linking in the preparation 
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nLo/^'^-nti/rf^-expl-i/j-f 
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dr 
L 
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(9.4) 



We observe that /gc is invariant under transformations of the auxiliary fields that correspond to independent trans- 
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lations and rotations of the replicas, as in the case of the 
Landau- Wilson free energy for the dimers-and-springs 
model (4.2). In addition, flf^ is a HRS field, f7^ is a 
IRS field, and Q" is a traceless and symmetric second- 
rank tensor field, all as with the corresponding auxil- 
iary fields for the dimers-and-springs model. Further- 
more, terms developed by expanding the log-trace part in 
Eq. (9.3) are structurally equivalent to those arising from 
the corresponding expansion of the Landau- Wilson free 
energy (4.2) of the dimers-and-springs model. In other 
words, as fields theories. Model A and the dimers-and- 
springs model are identical, up to elementary rescalings of 
their coefficients. We thus conclude that at length-scales 
larger than the microscopic scales at which they are de- 
fined, these models yield predictions that are structurally 
identical. 



ergy for Model B, a model for main-chain nematic poly- 
mer networks. This model consists of jointed chains (see, 
e.g., Ref. [39]) comprising S* — 1 rods each of length 
i. In addition, adjacent pairs of rods interact via a 
bending energy that promotes their parallel alignment, 
and arbitrary pairs of rods also interact via orientational 
forces that favor parallel or anti-parallel alignment. The 
network is formed via the random, instantaneous cross- 
linking of pairs of jointed chains via permanent point- like 
cross-links, located at the ends of the rods. The chains 
are labeled by z (= 1, . • • , N) and the rod-rod junctions 
are labeled by s {— 1, . . . , S*); the position of rod end s 
on chain i is Ci^s, as shown in Fig. 6. The system is then 
specified by the following Hamiltonian: 



B. Model B: main-chain nematic polymer network 

We apply the procedure used in Sec. Ill and the pre- 
vious subsection to to derive the Landau- Wilson free en- 



„ N S-1 , N S 
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2^4 



(9.5) 



;j = l s,t=l 



The partition function corresponding to this system, sub- 
ject to M cross-linking constraints, is given by 



N S 



^MC « nn /c«c,,,e-^"°/^(5(|c,,,+i-c,^,|-£) 

M 
X{5(c^.,s^-C^'^..s'^■ (9-6) 



i=ls=l • 

M 
X 

e=l 



Using the replica method and a Deam-Edwards type of 
distribution for the quenched randomness, together with 
a Hubbard-Stratonovich decoupling, we arrive at the fol- 
lowing Landau- Wilson free energy in terms of the auxil- 
iary fields {f]^,f^^,Q^}: 
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where A" 



A - ?72((5"'° - 1)/(NSV"''^), and the Boltzmann average 



/MC 



is defined to be 



/MC 



riLo nLi JdC^--- exp (2F ELo EL2 (c?+i 

nLo nf=i / dc^ exp (^ ELo Efr^ «+i - 



- • (c? 



<-i))nLonfr/>5(|<^i-<l-£) 
nLonL"i'^(|<+i-c?|-£) 



(9.8) 




FIG. 6; Model of a main-chain nematic polymer consisting of 
S— 1 rods of fixed length i. The chain is labeled i (= 1, . . . , A*') 
and the ends of the rods are labeled by s (where s — 1, . . . , S). 
Rod-end s on chain i is located at Ci,s. 



The comments made at the end of the previous sub- 
section concerning symmetries of the Landau- Wilson free 
energy under transformations of the auxiliary fields hold 
for Model B, too. 



X. CONCLUDING REMARKS 

The objective of this Paper has been to develop a 
microscopic approach to the liquid crystalline proper- 
ties of isotropic-genesis nematic elastomers (IGNEs), in 
which local nematic order — both in the preparation and 
the measurement ensembles — and random localization — 
induced by the presence of an elastomer network — are 
naturally incorporated. This development, which takes 
as its starting point a system of dimers that are perma- 
nently connected at random by Hookean springs, serves 
as the underpinning to the phenomenological approach 
presented in Ref. [22] by providing a systematic deriva- 
tion of the formulas on which the phenomenological de- 
scription is based. 

Specifically, by deriving an effective Hamiltonian of liq- 
uid crystallinity in IGNEs, we have shown that this mi- 
croscopic approach leads to the phenomenological Lan- 
dau theory of IGNEs proposed in Ref. [22] which, inter 
alia, predicts that at sufficiently large disorder strengths, 
both the thermal and glassy spatial correlations of ne- 
matic alignments can undergo oscillation with decay. 
The development has as a core feature an ensemble — 



the preparation ensemble — that is distinct from the usual 
(measurement) ensemble of Gibbs statistical mechanics. 
This feature enables us to determine in detail the influ- 
ence of the preparation history of IGNEs on their subse- 
quent equilibrium behavior. The consequences of these 
two ensembles were analyzed in detail in Ref. [22]. How- 
ever, the appearance of the two ensembles took the form 
of a hypothesis in Ref. [22], whereas in the present Pa- 
per they come into play naturally. Lastly, we have ar- 
gued that at sufficiently large length-scales, predictions 
made on the basis of a simple dimers-and-springs model 
are qualitatively identical to those made on the basis of 
two more realistic (but more complicated) microscopic 
models of IGNEs: one for side-chain liquid crystalline 
elastomers and one for main-chain liquid crystalline elas- 
tomers. 

Apart from its relevance to the specific subject of liquid 
crystalline elastomers, the present work brings to light a 
more general issue, viz., that the concept of a quenched 
random field should be broadened to incorporate not only 
the conventional, "frozen" type, which does not fiuctuate 
thermally, but also the type necessary for understand- 
ing media such as liquid crystalline elastomers, in which 
the frozen nature of the random field is present only at 
longer length-scales, fading out as the length-scale pro- 
gresses through a characteristic localization length, owing 
to the thermal position fluctuations of the network's con- 
stituents. The framework elucidated in the present work 
can be extended, with suitable modifications, to explore 
the statistical physics of other randomly cross-linked sys- 
tems, such as smectic elastomers and various biological 
materials. 
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Appendix A: Calculation of the average linking number per dimer 

In this series of Appendices we show various components of the calculations that are necessary for deriving the 
results presented in the main body of the Paper. In the present Appendix we show that the average number of springs 
connected to a dimer, [M]/7V, is approximately given by 2jf . By definition, 



[M] = E^(X)M 

X 

= J2 fvc"Pix\{c^})P{{c^})M 
X •' 

oo N N 
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•^ „-,- i„+ 11 „■ 1 



2iV(27r62)C/27" , 



To find a rough estimate for [M]/iV, we shall approximate the system as a dilute gas of dimers, which involves 
approximating the quantity exp(— jc^ — Cj|2/(262)) by (27r62)^/2^y_ This leads to 



N ^ N 



,(27r&2)^/2 



^ — 1 ^ + — 11'' ^- — 1 



Substituting this into Eq. (Al), we obtain 



[M]/iV«2r?2, (A3) 



which shows that rf provides a quantitative measure of the average number of cross-linking springs per dimer. 
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Appendix B: Hubbard-Stratonovich scheme 

In this Appendix, we present the details of the Hubbard-Stratonovich scheme that leads to Eq. (4.2). We first 
consider the replica partition function Zi^n in Eq- (3.10), which can equivalently be expressed as follows: 



oo N N ^ / -,r 9 \ M 



Zl+n - Y. H ■■■ H H ■■■ H W\ 



M=Oil.Jl = l iM,JM = l Si,ti = l,~l SM-tM 
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— \ \ P 2b2 ^^a = I t,s j,tl — p 262 I i,s j.f I 



"" '^''P UiV(27r62)^/2 E E 
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n2 

n N 
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Note that the two terms in the first exponential factor in the first equality are replicated versions of the Maier-Saupe 
interaction -ffncm and Hookean spring potential -ffev Also note that the first term in the last exponential factor in 
the first equality, - (T/r/2/2iV(27r52)^/2) ^^^^. ^^^^ exp (-|cO, - 4,|V262), jg ^^^ ^^^^ i/norm({< J)/rO in Eq. (3.5), 
whilst the second term in the same exponential factor is a replicated version of the excluded- volume interaction, ifov 
In the last step we have summed over all realizations of the quenched disorder. This results in an exponentiation of 
the Hookean energy term, and leads to the following expression: 
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(B2) 



Here, we have replaced Jij by its continuum limit J{cf — c"), viz., Ji,j « (y/A^) J(cf — c"), where we approximate 
J(c) by a potential of Gaussian form, i.e., J{c) « (Jo/(27ra2)''/2) exp(— 0^/20^), suitable for describing short-range 
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interactions. Next, we define the following collective fields in Fourier space: 

N 



qldM ^ ^Ee^^"-°Kd,<d.-^"''5d,,J, (B3a) 

1 ^ 

i=l s=-l,l 

1 ^ 

-^ - ^E E ^^""'■^- (B3c) 

i=l s=-l,l 

In Eq. (B3c), the argument of w^ can take any value in replicated Fourier space. However, in Eq. (B5) below, we 
shall decompose the field w^ into the part whose argument takes values from the HRS, the part whose argument takes 
values from the IRS, and the part whose argument takes values from the ORS. We shall then redefine Wj, to be only 
those parts whose arguments belong to the HRS. Note that Wpgc = w" and that the ORS part of ojj. is a constant. To 
simplify the notation, we define the following normalized averages: 

!NA+n = 



^. — n ^ — 1 ^ 



a={)i=\ 






47rt/^2 

We substitute the collective fields into Eq. (B2), taking care to separate out the ORS, IRS and HRS parts of the terms 
involving cj^,. The replica partition function then becomes 

z,,„ ^ iA^verl exp (^ E E ^i<..(p)P - f^ E ^1-0 P - ^ E E'^" 

^ \ Q— p a— 

X exp f -^Ic^ol' + ^ E E ^p\^P^- I 
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(B5) 



-p|yEEi&K^.(p)l'-yEE'^K"l 

Q=0 p Q=0 p ^"^ fe I ' N,l+n 





where ^ denotes a sum over all wave-vectors p excluding p = 0, we have made the following definitions: A^. = 
exp(-i62fc2)^ fji = 4^2^ ^^^ A^/T" = (4iV/l/)A"/r"-(4?7VT/")(l-(5"'0)Ap, and A^ is the effective excluded-volume 
interaction between dimers. For a 7^ 0, the excluded-volume interaction is renormalized downwards; owing to the 
presence of attractive interactions induced by cross-linking, whereas for a = 0, the excluded-volume interaction has 
no such renormalization. In the case of replica a — Q, the downward renormalization due to cross-linking is canceled 
exactly by the correction from the term -ffnorm! physically, this is expected as the dimers cannot feel the cross-link 
induced attractive forces at the instant just prior to cross-linking. 

We now implement the Hubbard-Stratonovich transformation, which is based on the following set of equalities for 
complex variables q and lo: 

e-^l«l' = - /"rf(Rew)d(Imw)e--^l'^l'+2'^^°«'^*; (B6a) 

e+-^l«l' - - /"d(Rew)d(Imc^)e-^l"l'+2-™^9-*. (B6b) 

Here, q may be thought of as a complex- variable analogue of a collective field and ui is the complex- variable analogue of 
its conjugate auxiliary field. The Hubbard-Stratonovich procedure allows us to express the replica partition function 
Zi^n as a functional integral in terms of the auxiliary fields fJ [whose argument we restrict to the HRS via the 
constraint l^peo, in accordance with the replica-sector division laid out in Eq. (B5)], 17" (which is formally il but with 
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its argument taking values in the IRS), and Q" [conjugate to (resp.) co, oj", and q"]. In terms of these auxihary 
fields, the replica partition function Zi+„ is given by 



/n 
vn II vn-^vct cxp ( - Nh+^ [n, n", q"]) . (b?) 

a=0 



Here, the Landau- Wilson free energy per dimer /i+„ (scaled in units of T) is given by 



^ a— p a— p 
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Q — P S— 1,— 1 
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By expanding /i+„ for small values of the auxiliary fields, one obtains a Landau theory in terms of Jl, J7" and 
Q", which are, respectively, the order-parameter fields for the random solidification transition, the crystallization 
transition, and the isotropic-nematic transition. As we assume that IGNEs are incompressible, there will be no 
fluctuations in the density of dimers (which can be enforced by making A to be extremely large), and thus there will 
be no corresponding instability in the IRS. We shall therefore disregard the contribution from $7". The free energy, 
Eq. (B8), then becomes 



~2 1 " J 

/!+«("' Q) = ^E^^l%l' + 2 E T.^\Q^dM? ln(cxp (Gi(fi) + G2(Q)))l.l+n . (B9) 
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Note that we have introduced the abbreviations 



Gi(") ^ W^H^k^u E ^"^'■"' (BIO) 
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G,{Q) ^ E E Y^QldM^^'""^ ( «. ^S,,,, ) ' (Bll) 
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and specialized to three spatial dimensions (i.e., D — 3). 



21 
Appendix C: Proof that (fi^)/ = [(a;^)J and (Q^/ = [(q^)^ 

We now prove Eq. (3.13), by deriving a general expression for [{qp^) . . . {qpp)] [P <E {1, 2,3, ... }) in terms of the 
conjugate fields Qp and the Landau- Wilson free energy /i+n [cf. Eq. (4.2)]: 
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I / 

[(gpi)x---(gpp>x] = 
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where in the second step, we have multiplied numerator and denominator by n — P copies of the factor 

„ N N 

^X = / n ^'^^l ^'^^-l ^^P (- (^ncm({nj) + iJev({Q, J) + i/xli„k({c,,,s J)) /T'') [] ,5(|c,,i - C,,_i| - ^) . (C2) 

This produces a factor of Z" in the denominator, which goes to unity once the replica limit is taken. We now perform 
the average over realizations of quenched disorder: 
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(C3) 

where -B„ is a constant. Here, \a\, . . . , ap} C {1, . . . , n}, as we are interested in nematic correlators in the measure- 
ment ensemble. In the second step, we have summed over M to obtain an exponential function. To proceed further. 
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we note that the factors of g" can be generated by introducing, into the rephca partition function ^i+„, a source 
field S)p that is hnearly coupled to q"; we denote the resulting replica partition function by the symbol Zi+„[Jo]. 
By functionally differentiating Zi_|_„[Jo] with respect to io" one can recover the correlators of q" At the end of the 
computation, one has to take the source limit io — > 0. Following this procedure, one obtains 

[{qpj^ . . . {qp,)^\ - hm hm ^„^ e ^ ^^„^^ . . . ^^„^ e 

xe *= 

' Ar,l+n 

hm hm ^^ — ^e^""'"^^'^'.^^^' -' ( -——- . . . -^^^^e ^ ^p t" iiv i-vi 



xe 
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xe ^ ^ . (C4) 

/ N,l+n 

In performing this calculation, we have completed the square in q" generating extra quadratic terms 
27V j° j° {-^p -^-p}- The Hubbard-Stratonovich transformation can now be performed, with Qp — -^j-^p, ^p, and w^, 
being the auxiliary fields. Making use of the following relation (which one can prove) 
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as we were aiming to establish. On going from the second to the third line of Eq. (C6), we have multiplied numerator 

^liq 



and denominator by Zi^n, and we have also taken the limit n — > 0, which implies that Zi-|_„ — > Z^- [as one can see 



from Eq. (3.10) by taking the limit n -^ 0]. For the case P = 2, we recover Eq. (3.13). 
We can also derive Eq. (8.2) by using Eq. (CI) and taking the replica limit: 

\(a )---(a )] = lim J ^"Hg^o^" ^^ I^pi ^pp)^ 
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^ ((QpJ---Qp;)), (C7) 

where in the second step, we have approximated the functional integrals by the values that are stationary with respect 
to n and exclude density fiuctuations i7". 
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Appendix D: Terms in the Landau- Wilson expansion 

In this section we expand the log trace term in Eq. (B9) for sniaU fl and Q". This gives us 

k 

a— p 

-^{{Gi{nfG2{Q)^)l,l+n - {Gl{nf)l.l+n{G2{Qf)ui+n)- (Dl) 

1. Terms proportional to Qfl 

First, we compute the quadratic term for the vulcanization part of the Landau theory: 
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At length-scales much larger than I or h, we can approximate the above result by 
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2. Terms proportional to Q.Q.Q, 
We next compute the cubic term for the vulcanization part of the Landau theory: 
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At length-scales much larger than i and 6, we can approximate this result by 



3. Terms proportional to Q^Q'* 

Next, we compute the quadratic terms from the nematic part of the Landau theory: 

(G,(0)^)M+n = Y.Y. E ^2QX^MQ'k.M) 

a,l3 p,q s,t=l~l ^ ' 
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Note that this term vanishes for a ^ fi. For a = /3, we find that 



f 111 

Y^/./. iT^Y' ^'did.2iP)Q%dS^P){^dtd3<^d2di +SdidJd2d3 - 7^Sdid2Sd3di)- (D7) 

a p ^ -^ 



Here, the notation J dn denotes integration over the unit sphere, and we have used the cquahtics (vafid for D = 3): 

15 ' 



J dnndind2 = ^<5did2 and J dnndind2nd3nd^ =■ |f ((5did2<5d3d4 +<5did3^d2d4 + ^d^di,^d2d3)- 



4. Coupling terms 

a. Terms proportional to flQ^Q^ and QflQ" 

The terms that couple the ncmatic order parameter to the vulcanization order parameter arc the ones that give rise 
to the physics of ncmatic elastomers. At cubic order in ft and Q (i.e., ilQ^Q*^ or i7i7Q") there are two such terms. 
They were computed in Ref. [23] and are given by 
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^{pd,Pd2 + {k" + ^P)di(fc" + ^p)d2)- (D8b) 

b. Terms proportional to fiilQ^Q'^ 

Terms proportional to ililQ^Q'^ arise in two forms: connected and disconnected. We first consider the disconnected 
type, which is given by {Gi{fl)'^){G2{Q)'^)- By using the expressions obtained for (Gi(rj)^)i,i+„ in Eq. (D2) and 
(G2(Q)^)ia+„ in Eq. (D7) we obtain 

(Gi(0)2)i,i+„(G2(g)2)i,i+„ = ^EEE^l (^) ^k^-kQld2iP)Qd,d2i-p)- (D9) 

Next, we consider the connected type, (Gi(f7)^G2(Q)^)i,i+„: 

" 7 7 

{Gm'G,{Qf).s+n = v'Y: Efc,.,,E E ^kAk2^k.^k2^j^QidMQidAci) 

a,p=0 ' p,q s,t = l,-l 

^i« \^d,d2)inin^d, \^d3d,))^ ^^ ■ (DIO) 

o o / 1,1+n 

Note that in this expression we have to evaluate terms such as (e*P'^"')i.i+„, ('^di'^d2 c*^''^°)i,i+n and 
(^di 'T'd2 ^da '^d4 6*^'^=) i_i+„. Wc recall the definitions 

Cs = c+{s£/2)n, (Dlla) 
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Here, we have decomposed coordinates into center of mass c = (ci + c_i)/2 and relative u = Ci — c_i, and have 
defined the integral over the unit sphere as 

f^^lfjcosO, (D12) 

where 9 is an angle that we can define with respect to the direction of the wave-vector. The following formulas are 
essential for the next step in our calculation: 

dn^.p.„ ^ sin^ 
At: \p\ 



/ -^(i'^''^{nd^nd., - -5d^o 



dn f d d 1 



47r \dpdi dpd^ 3 

1 f _d d_ , 1 

2 \dpd^ dpd 



+ o^did^ e 



o«pn- 



+ \5d,d^ j d(cos0)e^lPl^°^^ 



9 9 1 ^ sin IpI 

+ o"did2 



dpdi dpd^ 3 'V IpI 
1 3 \ , , 3 , A / |p|2 



^(rw4^--)' (°^^^) 



^ in-n / sm|p| 3cos|p| 3 sm |p| \ . 

-nd^Ud^ndznd^e ^ — I j — p^ -—r 1 j — r= — ((JdidaOdgdi + 0did,0d2d4 + Odidi0d2d^ 

\ \p[^ |p|4 IpI'^ y 



dn^ „ Avn ( sin|p| 3cos|p| ^ 3sin|p 

47r 

cos IpI 6sin|p| 15cos|p| 15sin|p 

+ SdidiPd2Pd3 + Sd^d^PdiPdi + Sd^d^PdiPds + SdsdiPdiPd; 

sin IpI 10cos|p| 45sin|p| 105cos|p| 105 sin |p 



4 ' |„|5 ' |„|6 IjjJT J {Sdid2Pd3Pdi + Sdid3Pd2Pdi 



|p|5 |p|6 |p|7 

15 " 210 

105 1890 83160 

+ 5d2d3PdiPdi + 5d2diPdiPd3 + Sdsd^PdiPd 



)Cos|p| iUt)sm|p|\ 

-Jp\8 + ^^^^)Pd,Pd2PdsPd, 



+ ^F^ ) i^didiSdad^ + ^did3Sd2di + 5didi5d2d3) 

{Sdid2Pd3Pdi + Sdid3Pd2Pdi + 5dtdiPd2Pd3 



(A \p? , IpI" 

V945 20790 1081080 



jPd^Pd2Pd3Pdi- (D13c) 



26 



l+n- 



Using Eqs. (Dll) and (D13), we can now compute (Gi(ri)^G2((5)^)i. 

""V 56 + 8064 56 + 8064 J^d^d^dd^d, 

+ ( - Y^ + 1008 )^ (P + 9)^1 (P + 9)d3'^d.d4 

+^^ b + '?)<ii (p + g)d2 (p + q)d; {v + g)d4 

+ ^(P + 9 + 4fc")<i, (p + q + 4fc")d,(p + g + ^k'-)dAv + g + 4fc")d4) 

+ 225 -^ Z^Z^y^2^^fe^-fe-pe°-qe'5^fc^-fc-p£--qe'3'3dld2(P)'3<i3d4(^) 

X {klk2Xy,^ + liPd.kl + klp,,){qd,kl + klq^J + Ipd.Pd^qd.qd,) , (DM) 

in which the Q" matrices are constrained to be symmetric and traceless. 

Appendix E: Deriving the effective Hamiltonian for liquid crystallinity 

In this Appendix we derive the effective Hamiltonian for Uquid crystallinity, Eq. (6.1). The first step in the derivation 
is to set fl to its saddle-point value Cl in /i+„ of Eq. (Dl) whilst retaining fiuctuations of Q, and evaluate the terms 
(Gi(17)2G2(Q))i,i+„, (Gi(t7)G2(Q)2)i,i+„, (Gi(f])2G2(Q)2)ia+„ and (Gi(f2)2)i,i+„(G2(Q)2)i,i+„ with n ^ n. As 
we are considering length-scales large compared with b or i, we shall make the approximations Aj « 1 and Jp Ri Jq. 
(The wave- vector dependent parts of A^ and Jp would generate corrections of higher order in the wave- vector.) By 
further arguing that fluctuation modes with p > ^^^ can be neglected in fci^, Q) (see App. E5), we shall arrive at 
the effective Hamiltonian (6.1). 



This term is given by 



1. {Gl{nfG2{Q))l,l+n 



^4^2 "J 1 1 



(Gi(f7)2G2(Q))l,l+„ « ^EEE^%^-fe-pe-QS,d.(P)(P<i.P<i.+(fc" + -p)rf,(fc" + -p)d,). (El) 

a=0 p J 

The product ^^.^_j^_ go involves a product of two Kronecker deltas, 5j2^_ fc°.o^-p-x;"_ fc°.Oi that enforces macro- 
scopic translational invariance, and implies that p = 0. Equation (El) then becomes 

n -4 /;2 J 

(Gi(ii)2G2(Q))i.i+n « E E \^^k ^^kkX k2.- (E2) 



a=0 



5T 



k 



By decomposing the full wave- vector sum ^j, into three separate parts corresponding to contributions from the HRS, 
IRS and ORS, the right hand side of Eq. (E2) becomes proportional to 

EIE- E - E \^k^-kKM.- (E3) 



a=0 



k kelRS keORS/ 
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Here, the ORS part vanishes, and the IRS part also vanishes, because ilj^ is proportional to 5y;"_ fe^.o, which is zero 
if k belongs to the IRS. We are thus left with the full wave-vector sum; it is proportional to 



k k 



= 0, (E4) 

where in the last step we have taken the replica limit. On going from the first to the second step, we have changed to 
center-of-mass coordinate c and relative coordinate u, and on going from the second to the third step we have 
shifted the wave-vector fc''' — > fe"*" — m/(2^|) and integrated over k^ (7 ~ 0,1,..., n) and r. Thus, the term 
{Gi{VLfG2{Q))i,i+n vanishes. 



2. (Gi(n)G2(Q)')i^ 



l+n 



To lowest order in wave-vector, this term is given by 

a^fi p.q 

ift V- V- [ dz JlG 



«#/3 pq •' 



a^P p.q 



- 800 EE^^P'^iP'i^P^3Pd.e-''«-g2,,,(p)gl,,(-p) (E5) 

a^0 p 

3. (Gi(n)2)i,i+„(G2(Q)2)i.i+„ 
This term vanishes. To show this, it is sufficient to consider the value of ^ f2^, ^_j^- Using the formula (5.9) for Q^ 

k 

and the decomposition ^ = ^j, — X^feeii^s ^ X^feeoflS' ^^ have that 

k 
fc V fe kelRS keOBsJ 

2]_(1 + ^)t 






(27rnF 



^ln( .o_;n.D )g' + Q("')- (E6) 

After taking the replica limit n — > 0, the above result vanishes. Here, we have used the definition ^j, = , wi+njo / dk 

when performing the integral over k. In the last step, we have integrated over fc, r, and c. The term 
(Gi(rJ)2)i.i+„(G2(Q)^)i,i+n is therefore zero in the replica limit. 

4. {G,{nfG2{Q)^)i,i+n 

We now consider the contribution (Gi(f2)^ G2{Q)^)- Despite its being of higher order in fj, this term has long wave- 
length properties that are sufficiently hard, and it is therefore important for the purpose of studying the polydomain 
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state of IGNEs. It contains a part that is diagonal in replica space (i.e., proportional to Q" Q") and a part that is 
off-diagonal in replica space (i.e., proportional to Q" Q^ for a ^ (i). 

First, let us consider the replica-diagonal part. From Eq. (D14), with Vt set to its saddle-point value fj, and 
by making use of the fact that i^^i^_^_(„ , „v° '^ ^Yl ki,o^-Yl k'-p-q.Oi which implies that p + q = 0, we see 
that replica-diagonal terms with prefactors in powers of \p + q\ vanish. Thus the replica-diagonal contribution to 
{Gi{n)'^ G2(<5)^) (with n set to Ct) becomes proportional to the quantity. 



J2J2J2^k^-kQldMQldA~p)[{'^-j^'\k'' 



a—0 p 



2... .-.,2 , V|^a|4^ 



63 



Sdid3Sd2 



1 



-16( - ^ + ^^Ik'^n^klk^Md, + ^f'klklkgkl 



(E7) 



The term proportional to Sj^-^^^dd2d4, (without the wave-vector prefactors) vanishes, as we showed in Eq. (E6). On the 
other hand, the remaining terms should also vanish, as f2^f2_^ is isotropic in wave-vector space. The vanishing of 

the term in X]q=o X^fc^d ^d ^k^-k comes about by the same argument as was presented in App. El. By means of 
a calculation similar to that performed in App. El, one can also show that the following result holds in the replica 
limit: 



2^fcdifcd2^d3^d4^fc^-fc — Z^ Z^ Z^ P^dl^d2^d3^34^fc^-fe - 0- 



(E8) 



kelRS keORS/ 



Next, we compute the replica off-diagonal contribution to (G'i(fi)^ G2(Q)^), whose terms we divide into the following 
three classes: (i) a term with prefactor proportional to kkkk, (ii) terms with prefactors proportional to kk but not 
kkkk, and (iii) terms with prefactors that do not depend on k [see Eq. (D14)]. 
First, we compute the replica-off-diagonal term proportional to kkkk: 

n 

T.T.T.^k^-k-p^--<i^^QldMQidM)''l''dydA 

a#/3 P,q k 



■2-,rn+3 



&-V 



X exp 



<l 



Hi 



J-J- (27r)3 (27r)3 (27r)3 7 V V 2^ ^did^yPJ^d^d^H) '^d^Kd^i^d^'^di 






k'^ 



.Zl- Z2 



,a , .21-22 , 1 



.21-22 1 

^^er+2^ 



a 



(n-l)|2i- 


22 1^ 


21-22 1 

2ei +2^ 


2 



1 



mp +q)-i{p + Q)-z2 



kf' + i 



.21-22 , 1 



2c: 



(E9) 



We replace the coordinates Zi and Z2 by the relative coordinate R and centcr-of-mass coordinate C, defined respec- 
tively hy R = Zl — Z2 and C = {zi + 22)/2, and integrate over C, which results in a factor of {2TT)'^VSp^-q^o. Next, 
we integrate over q to enforce the equality q — —p. We then integrate out k'^ (for 7 7^ a,/3). The right-hand side of 
Eq. (E9) then becomes 



Q2Yn+2 /^\(n-l)-D/2 " 
2 ] Z^ 



(27r)("-i)^ V^L 
X exp f - -p^^l - 



a, 13=0 
( = 7^)3) 



dfe" dk^ dp r dR 
{2n)D (27r)^ {2tt)D J ^ 



^"Vi'^'' "(ifc°i'+ife'n^i 



QldMQidS-p) 



:(fc° 



2- -- 4^i 

!^_P\ (],c._^_P\ (uf3_^,P) (uf3_}^^P 
2 o/diV^ 9^2 o)d2y Of'i O'dA 0C2 



X 



2ef 



2'd2' 



n 



2 / d3 ' 



2er 



2ldi- 



(ElO) 



Because the exponent is invariant under rotations of i?, fc", k^ and p, any terms having prefactors proportional to 
RdiRd2j ^d ^d 7 ^d ^d or pdiPd2 will, respectively, become proportional to 5did2 on integrating over i?, fc", k^ and 
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p. However, on contracting with C^did^j these terms wiU yield no contributions, owing to the tracelessness of Q,did2 
Thus, the term (ElO) effectively becomes 



g2y«+2 / ^ 



" f ^k^^k^J^ fdR a 

(2^)(»-i)^ V^iy „4^o 7 (27r)^ (27r)^ (2^)^ i y ^d^d,KP}^d,dS P) 

X exp ( - ip^ei - ^" + ^1'^'' - (|fe"P + Ife^nei 

X ( ifi^4 {RdiRd3Pd2Pdi + RdiRdiPd2Pd3 + Rd2RdsPdiPdi + Rd2Rd4.PdiPd3J 

-' sRdlRd2Rd3Rd4 + Tl^PdiPdiPdsPdi) (Ell) 



16^2 ^. ... ... -.4 -^g 

Then, by integrating over R, k" and fe'^, we obtain the following result: 

n 

E T.T.^k^-k-p^''-,^^Qid2ip)QidSQ)kik2ydA 



'•i^=" pq h 



2^4 



" 1 

E E ('^dld3'^d2d4 + 5d,d3Pd2PdAl + ■^Pd,Pd2Pd3Pdiii)e'~^P^^^^QXd2iP) Qd^dS^P) ' (^^2) 



i Q,/3 = p 



By means of similar calculations, we obtain the following results: 



E EE%^-fe-P^-.«^Q2,..(p)QL,(g) {Pd,k'^d2+klPd2){qd3ki+klqd:) 

k 
n 

E E {'^5d.d3Pd2Pd,eL +Pd,Pd2Pd3PdAl)e~''''^^^'Qld2iP) QidS-P)^ (E13) 



•f^=° PQ k 



^i a, 13 = p 



E EE^fe^-fe-pe---qe'^^rfl<i2(P)<3ld4(9)?'dlP'i29d39d4 



-''=° p,q k 



G' E T.p^^P'i^Pd^P'i^'~'''^^^'QldMQidS-p)- (E14) 



,/3 = p 



Putting the pieces together, we then obtain 

(Gi(l^)^G.(Q)^>M+« - IS^EE^ 



2-4 n4: JL 72 2 2 



/^2 -4 /I4 " 72 

225f f 2^ 2^ T'^Tf 

^-L 0,^ = p 



i^^ 



X (<5did3<5d2d4 + 5Sd,dsPd2PdAl + -^PdtPd2Pd3Pdi)Q'i,d2{P)Qd3dS^P)- ^^^^) 



5. Effective Hamiltonian for liquid crystallinity 

The effective Hamiltonian for liquid crystallinity, i/i+„[{Q"}^^Q] in Eq. (6.1), can now be derived using the 
coefficients obtained in App. D. Defining iJi+„[{Q"}^^Q] by 

Hi+„[{Q"}S=o] = /Q(Q) + /c(aQ), (E16) 
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we have 



1 " 7 1 

a—0 p 

-^(Gi(0)2G2(Q)>i,i+„ - ^(Gi(^)G2(g)2)ia+„ 

P \ a=l / 

(n 1 " 1 

-T^EE^p{QpQV+2T E ^pIQpQ^p} , (E17) 

a— 1 p a,/3 = l / 

where Hp is defined as in Eq. (6.2). Eq. (E17) is Eq. (6.1). In performing this derivation wc have made two 
approximations. Our first approximation consists in neglecting terms of cubic or quartic order in Q, as we are 
considering an IGNE that is both prepared and measured at high temperatures. Our second approximation consists 
in neglecting those terms in Eq. (??) that vanish with p^^ whilst retaining those terms that do not vanish with p^L- 
It may appear that terms with factors of larger powers in p^^ would dominate over those with smaller powers of pS,L 
for large p^l! however, the same terms also come equipped with exponential-damping factors of cxp(— p^^£/2), which 
means that the large p^^ contributions can effectively be neglected. We are thus left with the small p^^ modes, and 
for these, the terms in fc{^, Q) that vanish with p^^ are obviously smaller than those that do not vanish with p^l- 
The effective Hamiltonian thus obtained allows us to predict a novel regime in which correlators C"^ and C^ undergo 
spatial decay and oscillation (sec Sec. VIII B). As we shall sec in the following subsection, the terms neglected by our 
second approximation would not result in a qualitative change to this prediction. 



6. Corrections arising from neglected terms 

In the previous subsection the effective Hamiltonian _ffi+„[{Q"}2^Q] was derived under the approximation that we 
neglect terms originating in /c(^, Q) that contain factors of p^^. In this subsection we consider the modifications to 
the thermal and glassy correlators C^ and C^ arising from these neglected terms. For simplicity, we consider the case 
of a system prepared at a high temperature, so that the nematic order parameter in the zeroth replica, Q*', may be 
set to zero. Taking the neglected terms into account, one has the effective Hamiltonian 



jYeff[{Q"}s.i] = ^EE(-^*+^p'){%"^p"}-^EE^"^''^^{^p^p} 

a p Q^/3 p 

a p a^l3 p 



2T 

a p 



^ E E (-1(1 + '-]-PrL)e-'^''''^ + -2/eie-''^«^)<^^ . 4, (E18) 



2^2 

a^l3 p 

where we have introduced i^i = {G^fj"^ J^/ 225) {£/^l)'^ and ^2 = (Gry^ Jg/800)(i?/^L)^. We have also made use of the 
fact that any symmetric, traceless tensor Qp can be decomposed (in a wave-vector-dependent manner) into its five 
independent modes (f>p, Xp and Wpi viz., 

Q^i'^.(P) = ^ ^^ + Pl^i. ^ ^'iM + Pld. ^ ^-^3^) + W^^aAp), (E19) 
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where the modes (f>p, Xp and Wp are given by 

, _ Pdi Pd2 



IP, |p|0*<i.(P). (E20a) 

Wdid2ip) ^ Pdids Pd2di QdsdAP) - ^{Pdr^d.^ Pdsdi Qdsdi + (t>p)Pdid2, (E20b) 

XdAp) = Pd^d2j^Qd2dAp), (E20c) 

and the projection operator Pj^^^ is defined via 

By using Eq. (E18) we can compute the thermal and glassy correlators for the five independent modes in Fourier 
space. The thermal fluctuation correlators are given by 

2T 

[{{WpW.p)~{Wp){W^p)}] = .^^^ . ^ , . ... , (E22a) 

At + Lp'' + I'll ^e 2P iL 

2T 

[{Xp-X^p)-{Xp)-{X^p)] = ,, , , , , ^_,,^ , ^^ ,^,, _,,,. , (E22b) 



At + Cp^ + viT-^ (1 + 10p2^ 



e 2; 



[(0p^_p)-(</)p)(</)_p)] = ^- /^ ^ — Y- (E22c) 

We note that a given mode's correlator undergoes oscillatory decay in real space if the least stable mode has a nonzero 
wave- vector. Equivalently, if we Taylor expand the denominator of that correlator in powers of the wave- vector about 
zero and find that the effective stiffness (i.e., the coefficient of the quadratic term) is negative then the correlator 
undergoes oscillatory decay in real space. We can check that the effective stiffness does indeed become negative at 
sufficiently large values of the disorder strength (which is proportional to vi) for the modes Wp and (pp, although it 
does not for Xp. The full thermal fluctuation correlator C'^{r) is the sum of the inverse Fourier transforms of the 
correlators for the modes in Eq. (E22). Thus, at sufficiently large disorder strength, C^ {r) will undergo oscillatory 
decay. 

To determine whether the glassy correlator C'-^{r) undergoes oscillatory decay at sufficiently large values of the 
disorder strength we apply a similar analysis, computing the glassy correlator for each mode: 

[At + £p2 + j,^T-^{l + 10p2^2)g-ip2C| 

[(0p)(0-p)J = -2- (E23c) 



At + £p2 + r-i(i/i(l + f p4^i)e-5P'«i + V2pHie- 



p^a 



The full thermal fluctuation correlator C {r) is the sum of the inverse Fourier transforms of the correlators for the 
modes in Eq. (E23). As for the thermal correlator, the denominators of the glassy correlators for the modes Wp and 
<j)p (but not Xp) can acquire minima at non-zero wave-vectors for sufficiently large disorder strengths. In real space, 
such nonzero- wave- vector minima imply that the correlator C'-^{r) undergoes oscillation as it decays. 

Appendix F: A conventional random field approach 

In this Appendix we provide the steps leading to the thermal and glassy correlators that one would obtain from a 
conventional random field approach. In this approach, an IGNE subject to a given realization of quenched disorder 
H can be described by the Landau free energy 

Hri^\Y,{At + Cp'){ClpCl-p}-Y,{HpQ^p}. (Fl) 

p p 
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Here, H is a quenched random field having a Gaussian distribution defined by mean and variance 

[{Hdid-iip) Hd^dii—p)}] = ^{Sdxd2Sd3di + Sdid2Sd3d4 + Sdid2^d3d4)- 

By using the repHca technique to eliminate the quenched disorder, one obtains the efi'ective Hamiltonian 

-i n _. n 

Q = l p a, 13=1 p 

The inverse of the Hessian matrix corresponding to H^ is given by (for a, /? = 1, . . . , n) 



(F2a) 
(F2b) 



(F3) 






cap 



A 



I a/3 



{At + Cp^ 



(F4) 



where repeated Cartesian indices di, d2 mean that they are to be summed over, and 1"^ is an n x n matrix with each 
entry having the value of unity. Equation (F4) yields for the thermal and glassy correlators the results 



5T 



^p 

^p - 


At + £p2 ' 
5TA 


(At + Cp-^f 



(F5a) 
(F5b) 



Note that the denominators of C^ and C^ , calculated via a conventional random field approach developed in this 
appendix, do not feature the length-scale-dependent function Hp that appears in the denominators of (and plays an 
essential role in determining the behavior of) the correlators (8.3) derived from the microscopic dimer-and-springs 
model. 
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